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ABSTRACT 



Inexpensive Global Location and Tracking Systems are currently being designed to 
provide the civilian market low-cost radio position determination. This paper discusses 
two possible designs. The first design employs 3 or 4 satellites, depending on whether 
altitude is known a priori, each transmitting continuous ranging signals. The user 
transceiver receives the ranging signals, measures the time differentials of the receipt of 
the signals and transfers this information to a control station via a satellite link. The 
control station computes the user position from this data and sends the position coor- 
dinates back to the user via another satellite link. In the second design, each user 
transceiver transmits a unique code to the control station via the 3 or 4 satellite links, 
again depending on whether the altitude is known a priori. The control station measures 
the time differentials of the receipt of the signals and determines the user position. This 
position information is then transmitted back to the user via a satellite link. 
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I. INTRODUCTION 



Today, the United States Department of Defense (DOD) is utilizing the XAVSTAR 
Global Positioning System (GPS) for accurate position determination. GPS employs 18 
geosynchronous satellites within 6 different orbits with 3 satellites in each orbit. GPS 
provides highly accurate position determination to an unlimited number of users world- 
wide, although the position determination accuracy provided to the civilian community 
is intentionally degraded due to national defense reasons. The degraded mode still ena- 
bles the civilian community to attain accuracies of 100 meters. 

In addition to the high number of operational satellites that are required to support 
the 24-hour world-wide GPS system, the user equipment cost is also quite high due to 
its complexity. Even though the user equipment cost has decreased considerably since 
GPS was first introduced, the civilian community can not afford the user equipment. 
To lower the cost of the position determination system so that the civilian community 
can afford to install position determination equipment in one's vehicle, a simpler system 
has to be designed and the total number of satellites deployed has to be reduced. 

This thesis presents two low-cost concepts of vehicle location and tracking using 
geostationary satellites. The first concept is borrowed from the concept of geostationary 
satellite navigation systems presented in [Ref. 1]. The second concept is a new concept 
that is presented in this thesis. For the systems to be utilized by the civilian community, 
the systems have to be lower in cost than GPS, provide reasonable position accuracy, 
and be reliable. To cover the United States, each system requires four satellites to de- 
termine longitude, latitude, and altitude, or three satellites if altitude is not measured. 
Like the geostationary satellite navigation system, the location and tracking systems will 
also deploy three or four geostationary satellites parked in orbits on the equator at a 
latitude of zero degrees and altitude of approximately 35,786 kilometers. The difference 
between the proposed location and tracking system and the geostationary navigation 
system is that the location and tracking system user does not determine his position but 
relies on the central control station to determine it. The satellites continuously and 
synchronously transmit a satellite-unique ranging code to the user on earth. The user 
receives the satellites signals at different times due to the path length differences between 
the user and each of the satellites. The transceiver measures the time differentials be- 
tween a designated signal and the remaining signals and transmits this information to a 
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control station via a satellite communications link. The control station calculates the 
user position from the time differentials and transmits the user coordinates to the user 
via another satellite channel. By placing the computational burden on the central con- 
trol station, the location and tracking transceiver can be made at a lower cost. 

Position determination in a geostationary coordinate system is normally defined in 
a geocentric coordinate system as shown in Figure 1. The x-axis is the intersection of 
the equator plane (0° latitude) and the Greenwich meridian plane (O’ longitude) and is 
oriented from the center of the earth, the z-axis is the polar axis oriented from south to 
north, and the y-axis completes the right-handed Cartesian coordinate system. We des- 
ignate the unknown user position by (.v.y, z), the known positions of n = 3 or 4 
geostationary satellites by (x„y„ z,) where /= 1, 2, and the known position of the 
central control station by (.y 0 ,j’ 0 , z 0 ). 




Figure 1. Geocentric Coordinate System: [Ref. 1] 



2 



The main function of the transceiver is to measure the differential distance 



a i — d\ — d h i= 1,2,..., « (1) 

representing the differences in the distance to the user from one satellite, arbitrarily 
designated as satellite one, and from each of the other satellites. The variable a t repres- 
ents the difference in the distance of satellite one to the user and that of each of the other 
satellites to the user. The variable d, represents the distance from each satellite to the 
user. The differential distance determination concept is discussed later. The squared 
distance from the user position to each geostationary satellite can be expressed by 

d] = (x - xf + ( y -yf + (z - zf, i= 1,2, ..., « (2) 

where (x.y, z) represent the users coordinate position in the Cartesian coordinate system 
and (.Y rt y t . z) represents each of the satellites' positions. 

For n = 4. (1) and (2) combine to form a set of four independent equations and four 
unknowns, namely, x, y, z, and d x . These four equations are 

(4-a0 2 -(*-jr/+(y-.>v) J +(*-z i ) J , .'=1,2,... ,4 (3) 

which can be solved using the measured data a, and the known satellite ephemeris 
pt„y„ z), i = 1. 2. 3, 4. The value a, is found through the use of a correlation register and 
ranging codes which will be discussed later. A detailed discussion of the concept can be 
found in [Ref. 1], 

For the three satellite system, n = 3, we use 

x 2 +y 2 + ? = (R + hf (4) 

where R = 6378 km is the Earth's mean equatorial radius, and h is the user altitude 
(assumed to be known), as the fourth equation in addition to the three equations 
(4 - a) 2 = (x- x,) 2 + (y - y,) 2 4 - (z - z,) 2 i = 1 , 2, 3. 

The first system requires each satellite to broadcast a unique pseudo-noise (PN) 
code in the frequency band 2483.5 - 2500 MHz. The PX codes are selected from a set 
of Gold codes which allow minimum interference between signals from various satellites. 
The satellites employ extremely stable on-board atomic clocks which are mutually syn- 
chronized. A central control station which can view all satellites is responsible for the 
clock correction via a communications link in the frequency band 6526.5 - 6541.5 MHz. 
The user transceiver measures the time differentials between one satellite arbitrarily 
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designated as satellite one, and each of the other satellites. The accuracy of the meas- 
urements and their effect on the user position determination accuracy will be discussed 
later. The user transceiver then sends the measured time differentials to the central 
control station for computing the user coordinates via a code division multiple access 
(CDMA) uplink in the frequency band 1610 - 1626.5 MHz. The central control station 
calculates the user coordinates using the received time differentials and the predicted 
satellite ephemerides. It then relays the coordinates to the user via a time division mul- 
tiplexed (TDM) downlink in the frequency band 2483.5 - 2500 MHz. 

The second system requires the user transceiver to broadcast a PN code to the cen- 
tral control station via a code division multiple access uplink. The central control sta- 
tion calculates the user differential distances to the satellites with the available satellite 
ephemeris and then the user coordinates which are relayed to the user via the TDM 
downlink. 
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II. SYSTEM I 



A. SYSTEM I DESCRIPTION 

The frequency assignment and bandwidth allocated for each of the satellite links 
listed below are assigned by the Federal Communications Commission as listed in [Ref. 
2]- 

• Satellite-to-User Ranging Link 

Frequency 2483.5 - 2493.73 MHz 
Bandwidth 10.23 MHz 

• User-to-Control Station 

User-to-Satellite Uplink 
Frequency 1610 - 1626.5 MHz 
Bandwidth 16.5 MHz 
Satellite-to-Control Station Downlink 
Frequency 5150 - 5166.5 MHz 
Bandwidth 16.5 MHz 

• Control Station-to-User 

Control Station-to-Satellite Uplink 
Frequency 6524.5 - 6530.77 MHz 
Bandwidth 6.27 MHz 
Satellite-to- User Downlink 
Frequencv 2493.73 - 2500 MHz 
Bandwidth 6.27 MHz 

• Control Station-to-Satellite Command Link 

Frequency 6526.5 - 6541.5 MHz 
Bandwidth 15 MHz 

B. SATELLITE-TO-USER DOWNLINK 
1. Downlink Frequency and Bandwidth 

The communications downlink from each of the satellites to the user occupies 
a bandwidth of 16.5 MHz from 2483.5 to 2500 MHz. Within this bandwidth, the satel- 
lites must transmit a pseudo-noise ranging code needed for user position determination 
and provide the downlink for the control station-to-user data link. In System I design, 
the satellites must employ two satellite-to-user channels, one for the sole purpose of 
transmitting the unique pseudo-noise code and one to perform as a relay for the data 
link from the control station to the user. The pseudo-noise code will be transmitted 
between the frequencies of 2483.5 and 2493.73 MHz in a bandwidth of 10.23 MHz. The 
control station-to-user data link will be transmitted between the frequencies of 2493.73 
and 2500 MHz in a bandwidth of 6.27 MHz. 
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2. Pseudo-Noise (PN) Ranging Code 

The PN ranging code transmitted by the satellites is received by the user and is 
used to determine time differentials. The accuracy of the measurements of the time dif- 
ferentials used to determine the user position is of prime importance to the overall po- 
sition accuracy. The more precise the time measurement is, the more accurate the 
position determination will be. GPS uses code chip rates of 10.23 Mbit/s in its trans- 
missions to the user. For our example, a given code length of 2 19 — 1 == 524,287 chips 
and a chip rate of 10.23 Mbit/s results in a code period of 5 24,287/10.23 x 10 6 = 51.25 
ms. A lower code length of 2 16 —1 = 65,535 and a lower chip rate of 1.023 Mbit's has a 
period of 65,535. 1.023 x 10 6 = 64.06 ms and a larger error in the differential distance 
measurement. The code tracking loop that would be used is able to maintain code 
alignment to approximately 0.01 chip. This is a conservative figure, for code tracking 
loops can be accurate to 0.001 chip. For a 0.01 alignment factor and a 10.23 Mbit/s chip 
rate, a ranging error of 0.29 meter would exist, where a 1.023 Mbit, s chip rate would 
have a ranging error of 2.9 meters [Ref. 1]. To better understand the fundamentals of 
PN ranging codes that determine position accuracy, Figure 2 shows the concept of dif- 
ferential distance measurement. To determine the differential distance, the user equip- 
ment must be able to perform a correlation operation among the satellites' codes and 
internally stored copies of those codes. The receiver clock must be able to count the 
number of chips between the correlation peak of the first received code and the subse- 
quent codes as shown in Figure 2. Since a, may not be an integer multiple number of 
chips, the receiver may use the correlation value relative to the peak correlation voltage 
to determine a fraction of a chip in the measurement. As was stated earlier, a fraction 
of a chip of 0.001 can be measured. 



6 




Figure 2. Differential Distance Measurement: [Ref. 1] 

In [Ref. 1], the PN ranging code and its effects on position determination accuracy is 
discussed in further detail. 

C. USER-TO-CONTROL STATION LINK 
1. Uplink Bandwidth 

The frequency bandwidth allocated to the user-to-control station data link is 
1610-1626.5 MHz. The user direct sequence spreads the data with a Gold code of period 
2 10 - 1 = 1023 chips which is a standard period for navigation systems. Thus for a 
bandwidth of 16.5 MHz from 1610-1626.5 MHz, the typical user data rate is approxi- 
mately 16,129 bits.'s when spread at 1023 chips per bit, if Quaternary Phase Shift Keying 
(QPSK) is used as a modulation scheme. That is, the chip rate is ( 16,1 29)( 1,023) = 

16.5 x 10* chips. With QPSK, there are 2 bits per symbol and hence the symbol rate is 

8.5 x 10* symbols a second. The bandwidth required (double sideband for QPSK) is 
twice the symbol rate or equal to 16.5 MHz which is compatible with the available 16.5 
MHz bandwidth. The users must share this available bandwidth with each other to 
transmit their IDs, time differentials, and other information to the control station. 

The user-to-satellite uplink is the most critical link in the system. The user is 
not able to employ or utilize a high gain dish antenna like those used in the fixed 
stations. The maximum output RF power is also limited by the power available in a 
mobile vehicle. To calculate the Effective Isotropic Radiated Power (E1RP) of a vehicle, 
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Eq. (5) is used, where P, is the transmitter power out and O t is the effective gain of the 
antenna. 



EIRP = P t G, (5) 

In our design, a transmitter power out of 80 watts which is equal to 19 dBW and an 
omni directional antenna with a gain of 4 dB will yield an EIRP of 23 dBW. These design 
values are selected with some judgment and are also being used in current navigation 
systems. 

The Free Space Path Loss (FSPL) is calculated using Eq. (6). 

FSPLidB) = 10 log( ) 2 (6) 

where /is the carrier frequency, d is the slant range, and c is the speed of light. At a slant 
range of 35,786 km and a carrier frequency of 1618 MHz, a FSPL of approximately 190 
dB is obtained. 

A standard satellite G,'T of + 3 dB K has been chosen for this design model and 
Boltzmann's constant k is given as -228.6 dBW/K-Hz. To calculate the uplink carrier- 
to-noise density ratio, Eq. (7) is used. 

-y (dB) = EIRP(dBlV) - FSPL(dB) + y (db,'K) - k(dBW:K-Hz) (7) 

The EIRP of 23 dBW, FSPL of 190 dB, satellite G/T of +3 dB/K and Boltzmann's 
constant of -228.6 dBW K-Hz are entered into the Eq. (7), obtaining an uplink carrier- 
to-noise density ratio of 64.6 dB-Hz. This density ratio is sufficient to support the data 
rate that the system requires for successful operation. 

For the downlink, the same calculations are performed to obtain the downlink 
carrier-to-noise density ratio. A FSPL of 198 dB is calculated for the carrier frequency 
of 5150 MHz. A satellite EIRP of 36 dBW and a control station G'T of 20 dB K are 
used to calculate the downlink carrier-to-noise density ratio. These values are standard 
satellite EIRP and station G/T values. 

The total carrier-to-noise density ratio is calculated using Eq. (8) found in [Ref. 
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To solve for the total carrier-to-noise density ratio, the values of the carrier-to-noise 
density ratios must be converted to real numbers, substituted into Eq. (8), the equation 
solved and the answer converted to the logarithm form of dB-Hz. 

A margin of 2 dB is further engineered into the system design to allow for in- 
terferences that would degrade the operation of the system. The data rate is calculated 
by subtracting the 2 dB margin from the total carrier-to-noise ratio, converting that dB 
value to a factor and then dividing that factor by the equivalence of the required 
signal-to-noise ratio which is set to two representative values of 15 and 20 dB for this 
system design. These required signal-to-noise ratios are typical satellite navigation 
signal-to-noise ratios; each will be used for design calculations. 

To calculate the data rate for a lower required signal-to-noise ratio of 15 dB, the 
2 dB margin is subtracted from the total carrier-to-noise ratio of 64.57 dB-Hz, giving 
62.57 dB-Hz. This ratio is converted to a value of 1.807 x 10 6 and is divided by 31.26 
(value of 15 dB) which in turn yields the data rate of approximately 57 kbps. When the 
required signal-to-noise ratio is increased to 20 dB (factor of 100), the data rate is low- 
ered to IS kbps. The link analysis is tabulated in Table 1. 

2. User-to-Control Station Data 

After receiving 3 or 4 satellites' ranging signals and measuring time differentials, 
the user must transmit the information to the control station for calculation. The 
transmission is Code Division Multiple Access (CDMA) as required by the Federal 
Communications Commission (FCC) and is designed to provide the system with a bit 
error probability of 10~ 6 to 10~". (A bit error probability of 1(F 6 and 10 - ' provides the 
reliability needed for a successful system.) The data stream must include a synchronous 
code word in the user transmission for the control station to synchronize to the trans- 
mission. Once the the control station has synchronized itself to the data transmission, 
the control station can identify the information that the user is sending to the control 
station from within the data frame. A Barker code word of 13 bits would provide the 
synchronization that the control station equipment requires and enable the control sta- 
tion to lock on to the transmission and recover the user transmitted data [Ref. 4). 

Other information that is sent includes the user ID, the ID of the satellite that 
was received first, the second and third satellites IDs and the time differentials AT in 
symbols between receipt of the signals of the satellites. These time differentials are de- 
termined by the correlator. The minimum distance that a satellite may be from the user 
is 35,786 km which corresponds to a minimum propagation delay of 119.287 ms. The 
maximum distance is 41,678.82 km which corresponds to the maximum propagation 



9 



Table 1. SYSTEM I USER-TO-CONTROL STATION LINK ANALYSIS 



Uplink Frequency 


1618 MHz 


Antenna Gain 


4 dB 


Transmitted Power 


80 W 


Carrier EIRP 


23 dBW 


Free Space Path Loss 


190 dB 


Satellite G T 


+ 3 dB K 


Boltzmann's Constant 


-228.6 dBW K-Hz 


Uplink 

Carrier-to-Noise 
Density Ratio 


64.6 dB-Hz 


Downlink Frequency 


5150 MHz 


Satellite EIRP 


36 dBW 


Free Space Path Loss 


198 dB 


Control Station G T 


20 dB K 


Downlink 
Carrier-to-Noise 
Density Ratio 


86.6 dB-I Iz 


Total 

Carrier-to-Noise 
Density Ratio 


64.57 dB-Hz 


Margin 


2 dB 


£ 

If required = 


15 dB 


20 dB 


Then Data Rate can be 


57 kbps 


18 kbps 



delay of 138.929 ms. A maximum time differential between the minimum and maximum 
time delays is 19.642 ms. We require an accuracy of 0.1 n sec. This allows a quantization 
interval of 0.2 /isec. The number of 0.2 /xsec intervals in 19,642 /xsec is 98,210. To rep- 
resent this value, 2 18 = 262,144 will easily suffice to represent the maximum time differ- 
ential. Thus there are 18 bits required to represent the maximum differential time 
measurement between the first satellite and each of the other satellites and 2 additional 
bits to identify the satellite for a total of 20 bits used. In the three-satellite design, the 
altitude of the user must also be transmitted to the control station for position determi- 
nation. An altitude of 100,000 meters with 0.1 meter accuracy can be represented with 
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20 bits. For the four-satellite case, the field for altitude is used for satellite 4 ID and 
differential time measurement. In all cases, 20 bits for additional information are de- 
signed into the frame for future uses such as emergency reporting of position location 
when the user may be in need of it. See Figure 2 for an example of a data frame. 



ill BITS 



SAT 2 
ID & D T 



SAT 3 
ID & D T 



AL T I TUDE 



SAT 4 
ID & D T 



GENERAL 

INFORMATION 



Figure 3. Data Frame (User to Control Station) 



To increase the number of channels that could exist on the system, Forward 
Error Correction (FEC) is used. The use of FEC lowers the required signal-to-noise 
ratio and thus allows more CDMA channels to be used. A convolution code with 
Viterbi decoding and soft decision will increase the total number of channels [Ref. 5J. 
There are various combinations of convolution code rates and constraint lengths that 
can be examined to provide the optimum data rate for the system. 

To calculate the number K+ 1 of CDMA channels that can successfully exist 
within the bandwidth at a bit error probability of I0 6 and 10"\ the following results 
from Refs. 5 and 6 were used. To calculate the number K+ 1 of CDMA channels that 
can be succesfully utilized simultaneously, Eq. (9) is used. 



K = 3 1\ 



Jo_ 

2E h 



2 E h 



(9) 



The chips per bit (N), the design bit energy-tomoise ratio 
to-noise plus user interference density ratio ( —E \ must be 



*)■ 



and the bit energy- 
mown to solve for K. For 



reliability, bit energy-to-noise ratios of 15 and 20 dB were selected for the system. The 
chips per bit (X) is 1023. The value for the bit energy-to-noise plus user interference 
density ratio is a function of other parameters as follows. 




(10) 



where <2 (jc) is defined as 



ew- 



J2n 




(H) 



The code rate (R) and the channel transition probability (p) must be known to calculate 
the energy-to-noise p’us user interference density ratio. A code rate (R) of 1,2 to 7,8 
was used in the calculations to determine the bit rate that will support the maximum 
number of CDMA channels. The probability (p) is a function of parameters that are 
determined by the code rate and system design. To calculate p, Eq. 12 is used [Ref. 6]. 



(12) 



For System I design, coded bit error probabilities (P byCBded ) of 10 -6 and lO -1 are 
used to insure a reliable system for the user. The number k of encoder input bits and 
the number n of encoder output bits establish the code rate (R), where /? = — . The 
variable B d used to calculate the channel transition probability (p) is a calculated value 
that is determined by the code rate, constraint length, and the free distance of the code. 
The variable d fr „ is the free distance of the convolutional code and is determined by the 
code rate and the constraint length. The values B d/m and d fi „ have been previously cal- 
culated in [Refs. 5,6] at a constraint length of 7 for the code rates of 1 2 to 7/8. 

To calculate the number of simultaneous CDMA channels that the system can 
support, the channel transition probability (p) must be found with Eq. (12). The number 
k of encoder input bits, B dfree , d /m and the coded bit error probability (P btCOded ) are entered 
into Eq. (12) to solve for the channel transition probability (p). Then the bit energy- 
to-noise plus user interference density ratio is substituted into Eq. (10) along 

with the given chips per bit (X) of 1023, and the design bit energy-to-noise density 
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ratio ^ -p- of 15 and 20 dB. In Tables 2 and 4, the number K + 1 of channels that can 
successfully be used with hard decision Viterbi coding at a bit error probability of 10~ 6 
and 10 - ' respectively are shown. For example, given the code rate (R) of 1,2, the en- 
coded bit k is 1, B df " is 36, and d /r „ is 10. For a P Kco<ttd of 10 6 a channel transition 
probability (p) of .00703 is calculated. With a p of .00703 and a code rate (R) of 1,2, a 
bit energy-to-noise plus user interference density ratio ^ -jp ^ of 5.87 is calculated. This 
value, the chips per bit (N) of 1023, and the design energy-to-noise density ratio of 20 
dB (100) yields a value K of 246. The number of channels K + 1 is then 247. An in- 
spection of the tables reveals that a code rate of 1/2 yields the largest number of simul- 
taneous channels that can be used successfully. When soft decision Viterbi coding is 
used, a decrease in the bit energy-to-noise plus user interference density ratio of 2 dB is 
realized. Thus soft decision Viterbi coding further increases the number of CDMA 
channels that can be used which is shown in Tables 3 and 5 for 10~ 6 and 10“' respectively. 
The different bit error probabilities of 10~ 6 and 10 -7 and soft and hard Viterbi decisions 
are used in the calculations to provide the reader an indication of four possible system 
capabilities. 

The major limitation in using lower code rates of 12 is that the bandwidth of 
the system that is available to the user is severely reduced. Although the number of 
CDMA simultaneous channels is greatly increased, the bandwidth of the system is de- 
creased. Further research is required in the area of the coding techniques to find the 
maximum number of CDMA channels that can be used with the highest coding rate to 
allow for the greatest user bandwidth. 

3. Satellite-to-Control Station Downlink 

The frequency bandwidth allocated for the satellite-to-control station downlink 
is 16.5 MHz from 5150 to 5165.5 MHz. The satellite serves as a relay for the user-to- 
control station data link. 

D. POSITION DETERMINATION 

After receipt of the time differentials, the central control station computes the co- 
ordinates of the user from the measured time differentials, and the known location of 
each of the satellites. 

The x, y, z position determination of the user is defined in a geocentric coordinate 
system as shown in Figure 1. The x axis is the intersection of the equatorial plane and 
the Greenwich Meridian plane and is oriented from the center of the earth. The y axis 



13 



Table 2. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 



(Hard Decision) at Coded Bit Error probability of 10~ 6 



Code 

Rate 

R 


1/2 


2/3 


3/4 


4/5 


5/6 


6/7 


7/8 


k 


1 


2 


3 


4 


5 


6 


7 




36 


3 


42 


12 


92 


5 


9 


dfree 


10 


6 


5 


4 


4 


3 


3 


CDMA Channels K 
for = 100 or 20 dB 


247 


237 


185 


172 


157 


147 


145 


CDMA Channels K 
for -^r- =32 or 15 dB 


213 


204 


152 


139 


124 


114 


112 



completes the cartesian coordinate system with the x axis. The z axis is the north to 
south axis connecting the two poles. 

Although actual distance is not measured by the user equipment, the distance can 
be calculated from the time it took for the other satellites signals to reach the user versus 
the first satellite signal that was received. The time differentials are converted to differ- 
ential distances by multiplying the time differentials by the speed of light (c); approxi- 
mately 3 x 10 ? m s. The accuracy of the synchronization of the satellites transmissions 
plays an important part in the users position determination. If the satellites do not 
transmit at the exact same time, the time differentials will not be correct and thus neither 
will the distance differentials. 

The distance between the unknown user position (x, y, z) and each of the satellites 
(x, z,) where /= 1, 2, .... n can be calculated from 

d] = {x- xf + (y -yif + {z~ '=1,2 n. (13) 

We further define the distance measurements 

a i = d l —d i /=!, 2, .... n (14) 
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Table 3. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 



(Soft Decision) at Coded Bit Error probability of 10 -6 



Code 

Rate 

R 


1/2 


2/3 


3/4 


4/5 


5/6 


6/7 


7/8 


k 


1 


2 


3 


4 


5 


6 


7 




36 


3 


42 


12 


92 


5 


9 


dfree 


10 


6 


5 


4 


4 


3 


3 


CDMA Channels K 
for = 100 or 20 dB 


400 


384 


302 


281 


258 


242 


238 


CDMA Channels K 
E„ 

for — = 32 or 15 dB 


366 


351 


269 


247 


225 


209 


206 



as the difference in distance from one satellite designated as satellite one to the user and 
the other satellites distance to the user. 

For a four satellite case, n = 4, equations (13) and (14) are combined to form four 
independent equations and four unknowns, x, y, z, and d v From the known satellite 
positions and the time differences which have been measured by the user, the four un- 
knowns can be solved. Sign ambiguities can be corrected as long as the user knows what 
quadrant he is in. 

The coordinates of the satellites (jr„ y„ z,) can be calculated from its longitude 
eastward of zero degrees by 



jc ( = r cos 6[ 


(15) 


St = r sin 6j 


(16) 


z ( = 0 


(17) 



The term z, will remain zero as the satellites are in an equatorial orbit with a latitude of 
zero degrees. 
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Table 4. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 



(Hard Decision) at Coded Bit Error probability of 10 - ’ 



Code 

Rate 

R 


1/2 


2/3 


3/4 


4/5 


5/6 


6/7 


7/8 


k 


1 


2 


3 


4 


5 


6 


7 




36 


3 


42 


12 


92 


5 


9 


dfree 


10 


6 


5 


4 


4 


3 


3 


CDMA Channels K 
E s 

for — = 100 or 20 dB 


215 


200 


159 


147 


135 


123 


122 


CD VIA Channels K 
£, 

for — ~ = 32 or 15 dB 


1S1 


167 


126 


112 


102 


89 


88 



The four non-linear independent equations are 

A - af = ( x - xf + (j y-yf + z 2 i = 1. 2, 3, 4 (18) 

Equation (18) can be reduced to 3 linear independent equations by substituting equation 
(14) into equation (13) and subtracting the (/ = 1) equation from the other 3 equations. 
The remaining three equations are 



{xt-xjx + t vi-y\)y ~ a i d \ =~-y- *-2, 3, 4. (19) 

The 3 equations can then be entered into matrix form to solve x, y, and d x as shown: 



X 




*2-*i yi~y\ ~ a 2 




«2 




l 




-1 


2 


y 


2 


*3-* i y'3 -y\ ~ a 2 




a 3 


A. 




<3 

1 

>1 

1 

1 

>? 




a\ 



z can be calculated by substituting equation (20) into equation (18) 
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Table 5. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 



(Soft Decision) at Coded Bit Error Probability of 10~ 7 



Code 

Rate 

R 


1/2 


2/3 


3/4 


4/5 


5/6 


6/7 


7/8 


k 


1 


2 


3 


4 


5 


6 


7 


K. 


36 


3 


42 


12 


92 


5 


9 


dfree 


10 


6 


5 


4 


4 


3 


3 


CDMA Channels K 

El 

for -yr- = 100 or 20 dB 


349 


326 


261 


239 


223 


203 


201 


CDMA Channels K 
for = 32 or 15 dB 


316 


292 


228 


206 


190 


170 


168 



; = [(</,-(./ - (x — Xj) 2 - (y-y, ) 2 ]i. (21) 

From equation (20) and (21), the user position (x, y, z) can be solved to obtain an initial 
approximation. The actual position must be solved iteratively or by utilizing a Kalman 
filter. From (x, y, z), the users altitude h, longitude and latitude can be derived. 

x = {R + h) cos 6 l cos 6/ (22) 

y = {R + h ) sin d L cos Q t (23) 

z = (R + h) sin d t (24) 

where 0 /t and h are solved by 

Si-tan-'f-i) (25) 

cos 6, 

9i= tan \z — ^- L ) (26) 
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(27) 



For a three satellite system, the position solution algorithm is much the same. The 
following equation substitutes for the fourth satellite equation 



= (R + h? 



(28) 



As in the four satellite case, these four equations must be solved iteratively to obtain 
the users position. In the three satellite case, the user must report his altitude because 
only the longitude and latitude of his position can be determined from three satellites. 
The four non-linear equations are 

{d\ - af = (x - xf + {y-yf + 2 1 1=1, 2, 3 (29) 

x 1 +y 2 + z 1 = {R + hf. (30) 



As in the four satellite case, the equation with index = 1 is subtracted from the 
other two equations. 



+ B,^ + y j 


(31) 


d l = a 2 x + + y 2 


(32) 



where 



y*-y i *i + Ji - *2 - >2 + a i 

?i = - 



2a, 



X\ + Tl - *3 - >3 + a 2 
2 a. 



When we solve for y in terms of x, we get 



(33) 

(34) 

(35) 



y = H\X + H 2 



(36) 



where 
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(37) 



Mi = 



yi ~ 

B 2 -Bi ’ 



M 2 = 



y i - yi 
b 2 -b, • 



The substitution of (36) into (31) produces 

d\ = n 2 x + jU 4 



(38) 



where 



ot,B 2 a 2 B 1 ot[B 2 a 2 B] 

^ 1 B 2 - B ] B 2 - B, 



(39) 



The term z 2 can be expressed as 

z 2 = d]-x 2 - v 2 + 2 X] x + 2y\y - x 2 - y 2 . (40) 

Substitution of (36), (38), and (40) into (28) produces 

ax 2 + bx + c = 0 (41) 

where 

a -A (J2) 

b = 2m 3 m 4 + 2y ] n ] + 2 X] (43) 

c = /4 + 2y lf i 2 ~ x] - y 2 - (R + h) 2 . (44) 

The value of x, y, and z can be calculated from (36), (40), and (41) subject to the 
limitation x < 6378 km provided the user knows what quadrant of the earth he is lo- 
cated in. 

E. POSITION ERROR 

Even though the satellite is geostationary', the satellite does move about its nominal 
position. The drift of the satellite can be calculated to less than 20 meters from its 

nominal position. The error that results in the user position determination can be cal- 

culated. The range error can be expressed as user position error (Ax, Ay, A z.) 
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cd: cd-. cd. 

Ad-. = Ax 4- -r 1 Ay + Az 

cx cy cz 



-Ax 



y-y'i 



Ay + 



~ z t 



/= 1, 2, 3, 4. 



Az 



(45) 



d; d { ' di 

Equation (45) can be transformed into matrix form as Av = u and with equation (14), 



x — x l 


y-y i 


z-z. 


4 


dx 


d\ 


X — x 2 


y-y’2 


z-z 2 


d x - a 2 


d\ ~ a 2 


d \ ~ «2 


X — Xj 


y-y% 


z-z 3 


d x -a 3 


d\ ~ a 3 


d\ - «3 


X - x 4 


y-y 4 


z -z 4 


d x - o 4 


d\ ~ ^4 


d\ — a 4 



(46) 



(47) 



The user position error can then be simulated by letting the range error be 
A d t = 0~* # j) 2 + (y-y’f + {z-z’fj 2 



[(x-xf + (y-yf + (z-z,.) 2 ]T 



f'= 1. 2, 3, 4 



(48) 



where 



(x’b y'i’ z ’i) = (-Y; + Ax;, }>i + Ay h Zi + Azt). (49) 

The variable Ad, is the difference between the distance from the user to the nominal 
satellite position and the distance from the user to the satellite position that the satellite 
may have drifted to. The variables (x'„y' t , z' t ) are the coordinates of the nominal satellite 
position with the drift distances ( A.v„ Ay, A z, ) added to the nominal satellite position 
coordinates. The error vector v is the linear least-squares solution of Av = u where 
v = (A r A)~' A T u and the mean squared error is e = (Av — u) T (Av — u). The values 
(A.v, Ay, Az) are calculated using Eq. (46) and (47). The error contribution of the z co- 
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ordinate is much less than the x or y coordinate as the user position approaches the 
equator. This ill-conditioned problem results in a poor estimation of A z which can be 
minimized by modifying the least-squares method of solution. By first calculating the 
least-squares estimate of Ax and Ay and ignoring the A z component of the total error 
and calculating A 2 after Ax and Ay are found, a better estimate of the total error can be 
found. 



where 



». = [£] = (A[A,r'A[u 
A z = (A[a 2 )~ ] (A[u - a[a iV] ) 




(50) 

(51) 



(52) 



(53) 



For a three satellite system, the range errors Ad, i= 1, 2, 3 are calculated from 
equation (45) and the user altitude error Ah is calculated from equation (30) as 
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Ah = 



x(Ax) + j-(Aj-) + z(Az) 

[x 2 + > + z 2 ] T 
■x(A.y) + j(A>) + z(Az) 
R + h 



(54) 



The user position error vector v = CAx, A v v, A zJ T is the least-squares solution of 
Bv = w which is v = ( B T B )~ 1 B T w , where 



* - *1 y - y i z - z\ 



di 


d x 




d x 




x — x 2 


y - yi 


2 


~ z 2 




'a d~ 


d\ — a 2 


d \ ~ a 2 


d x 


~ <*2 




Ad 2 










w — 




x - x 3 


y - >3 


2 


~ *3 




A d 3 


d\ - a i 


d \ ~ «3 


d x 


- *3 




Ah 


X 


y 




2 




R + h 


R + h 


R 


+ h 





Tor the three satellite case, the least-squares method is also modified to provide a better 
estimate of v. The variable v is calculated from v = [v,, Az] r where 

= (<£,)“' B?w (56) 




A z = (B[ B 2 )~ l (B[ w - B 2 B x v,) 



x — x x 


y - J 




2 ~ 2 X 


d x 


d x 




d x 


X - x 2 


y - >’2 




z - z 2 


d X - a 2 
X - - r 3 


* % 

i i 


B 2 = 


dx - a 2 
* ~ 2 3 


d \ - a 2 


d X - a 3 




dx - «3 


X 


y 




2 


R + h 


R + h 




R + h 



(57) 



(58) 
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Equations (50), (51), (56), and (57) calculate the upper limit of the error in the user 
position provided the satellites position errors are known in either the three of four sat- 
ellite case. The average error for various users positions has been compiled in [Ref. 1]. 
The average error in the users position increases as the user approaches the equator. 
The average error ranges from 14.2 to 61.5 meters depending on how close to the 
equator the user is. The closer the user is to the equator, the worse the error which in- 
dicates a weakness in the system. 

F. CONTROL STATION-TO-USER DATA LINK 
1. Uplink Bandwidth 

The frequency bandwidth allocated for the control station-to-satellite uplink is 
6524.5 to 6530.77 MHz. The 6.27 MHz bandwidth coincides with the available band- 
width of the satellite-to-user downlink part of the relay. 

After the control station has completed the calculations and the user position 
is determined, the position coordinates have to be sent to the user. A TDM channel is 
used to transport the user information to the user. To calculate the total carrier-to-noise 
density ratio of the control station-to-user link. Eqs. (5) - (S) are used as was discussed 
in the user-to-control link. 

The uplink frequency is 6533 MHz with an antenna gain (G,) of 44.5 dB. which 
is a standard antenna gain for a ground station antenna. The antenna gain is much 
higher than for a mobile vehicle as one is able to use a large dish. The power out ( P ,) 
is set at only 1 watt per channnel since there would be more than one channel simul- 
taneously transmitted by the control station. For example, if 200 CDMA channels were 
used, then 200 watts of power would be required for transmission of the channels. The 
lower power out is standard for a ground station. For this link design, the satellite G T 
is -1.5 dB K. The power out and the satellite G/T values are design parameters deter- 
mined by judgement and current systems. For the control station-to-user link, the total 
carrier-to-noise density ratio is 66.25 dB-Hz is calculated by Eq. (8). A carrier-to-noise 
density ratio of 66.25 dB-Hz is sufficient for a successful link. A 5 dB margin is engi- 
neered into the system to allow for the cross-interference degradation of the system by 
the multiple carriers that are transmitted by the control station. After subtracting the 
5 dB margin from the total carrier-to-noise density ratio, converting to a factor and di- 
viding the result by the required noise ratio factor, the data rate for the bit error proba- 
bilities of 10~ 6 and 10*’ are 42.2 kbps and 13.3 kbps, respectively. The calculations for 
the control station-to-user link utilize Eqs. (5) - (8), which have been discussed earlier. 
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